Development of a Fully Automated Image Analysis Method for High Density cDNA and array CGH Microarray Based Genomic Studies

In this thesis, different techniques for image analysis of high density microarrays have been investigated. Most of the existing image analysis techniques require prior knowledge of image specific parameters and direct user intervention for microarray image quantification. The objective of this research work was to develop of a fully automated image analysis method capable of accurately quantifying the intensity information from high density microarrays images. The method should be robust against noise and contaminations that commonly occur in different stages of microarray development.

Caracterização de rearranjos cromossômicos citogeneticamente equilibrados associados a quadros clínicos

Este estudo teve como objetivos (a) identificar mecanismos pelos quais rearranjos cromossômicos citogeneticamente equilibrados possam estar associados de maneira causal a determinados quadros clínicos e (b) contribuir para a compreensão dos mecanismos de formação desses rearranjos. Para isso, foram estudados 45 rearranjos cromossômicos citogeneticamente equilibrados (29 translocações, 10 inversões e seis rearranjos complexos), detectados em pacientes que apresentavam malformações congênitas, comprometimento do desenvolvimento neuropsicomotor ou déficit intelectual. Foram 31 rearranjos cromossômicos esporádicos, três familiais que segregavam com o quadro clínico e mais 11 rearranjos cromossômicos herdados de genitores fenotipicamente normais. Inicialmente os pontos de quebra desses rearranjos foram mapeados por hibridação in situ fluorescente (FISH). A busca por microdeleções e duplicações genômicas foi realizada por a-CGH. A investigação dos pontos de quebra prosseguiu com a aplicação da técnica de Mate-Pair Sequencing (MPS), que permite localizar as quebras em segmentos de 100 pb - 1 kb, na maioria dos casos. Para obter os segmentos de junção das quebras no nível de pares de bases, os segmentos delimitados por MPS foram sequenciados pelo método de Sanger. A análise por aCGH revelou microdeleções ou microduplicações localizadas nos cromossomos rearranjados, em 12 dos 45 pacientes investigados (27%). A análise de 27 rearranjos por MPS permitiu a caracterização dos pontos de junção das quebras. MPS expandiu o número de pontos de quebra, detectados por análise do cariótipo ou aCGH, de 114 para 156 (em resolução < 2kb, na maioria dos casos). O número de pontos de quebra/rearranjo variou de 2 a 20. Os 156 pontos de quebra resultaram em 86 variantes estruturais equilibradas e outras 32 variantes não equilibradas. Perdas e ganhos de segmentos submiscroscópicos nos cromossomos rearranjados constituíram a principal causa ou, provavelmente, contribuíram para o quadro clínico de 12 dos 45 pacientes. Em cinco desses 12 rearranjos foram detectadas por MPS a interrupção de genes já relacionados à doença, ou provável alteração de sua região reguladora, contribundo para o quadro clínico. Em quatro dos 33 rearranjos não associados a perdas ou ganhos de segmentos, a análise por MPS revelou a interrupção de genes que já foram anteriormente relacionados a doenças, explicando-se, assim, as características clínicas dos portadores; outro rearranjo pode ter levando alteração da expressão gênica de gene sensível a dosagem e ao quadro clínico. Um rearranjo cromossômico familial, identificado na análise após bandamento G como uma translocação equilibrada, t(2;22)(p14;q12), segregava com quadro de atraso do desenvolvimento neuropsicomotor e dificuldade de aprendizado associados a dismorfismos. A combinação das análises por FISH, aCGH e MPS revelou que se tratava, na verdade, de rearranjo complexo entre os cromossomos 2, 5 e 22, incluindo 10 quebras. A segregação de diferentes desequilíbrios submicroscópicos em indivíduos afetados e clinicamente normais permitiu a compreensão da variabilidade clínica observada na família. Rearranjos equilibrados detectados em indivíduos afetados, mas herdados de genitores clinicamente normais, são, em geral, considerados como não tendo relação com o quadro clínico, apesar da possibilidade de desequilíbrios cromossômicos gerados por permuta desigual na meiose do genitor portador do rearranjo. Neste trabalho, a investigação de 11 desses rearranjos por aCGH não revelou perdas ou ganhos de segmentos nos cromossomos rearranjados. No entanto, a análise por aCGH da portadora de um desses rearranjos - inv(12)mat - revelou deleção de 8,7 Mb no cromossomo 8, como causa de seu fenótipo clínico. Essa deleção estava relacionada com outro rearranjo equilibrado também presente em sua mãe, independente da inversão. Para compreender os mecanismos de formação de rearranjos citogeneticamente equilibrados, investigamos os segmentos de junção no nível de pares de base. A análise por MPS que levou, na maioria dos casos, ao mapeamento dos pontos de quebras em segmentos <1kb permitiu o sequenciamento pelo método de Sanger de 51 segmentos de junções de 17 rearranjos. A ocorrência de blunt fusions ou inserções e deleções <10 pb, e a ausência de homologia ou a presença de micro homologia de 2 pb a 4 pb de extensão indicaram o mecanismo de junção de extremidades não homólogas (non-homologous end joinging; NHEJ), na maioria das 51 junções caracterizadas. As características de três dos quatro rearranjos mais complexos, com 17-20 quebras, indicaram sua formação pelo mecanismo de chromothripsis. Este estudo mostra a importância da análise genômica de variações de número de cópias por microarray, juntamente com o mapeamento dos pontos de quebra por MPS, para determinar a estrutura de rearranjos cromossômicos citogeneticamente equilibrados e seu impacto clínico. O mapeamento dos segmentos de junção por MPS, permitindo o sequenciamento pelo método de Sanger, foi essencial para a compreensão de mecanismos de formação desses rearranjos

Chromosomal gains and losses in ocular melanoma detected by comparative genomic hybridization in an Australian population-based study

To define the location of potential oncogenes and tumor suppressor genes in ocular melanoma we carried out comparative genomic hybridization (CGH) analysis on a population-based series of 25 formalin-fixed, paraffin-embedded primary tumors comprising 17 choroidal, 2 ciliary body, 4 iris, and 2 conjunctival melanomas. Twelve (48%) of the 25 melanomas showed no chromosomal changes and 13 (52%) had at least one chromosomal gain or loss. The mean number of CGH changes in all tumors was 3.3, with similar mean numbers of chromosomal gains (1.5) and losses (1.8). The highest number of chromosomal changes (i.e., nine) occurred in a conjunctival melanoma and included four changes not observed in tumors at any other ocular site (gains in 22q and 11p and losses in 6p and 17p). The most frequent gains in all primary ocular melanomas were on chromosome arm 8q (69%), 6p (31%) and 8p (23%) and the most frequent losses were on 6q (38%), 10q (23%), and 16q (23%). The most common pairing was gain in 8p and gain in 8q, implying a whole chromosome copy number increase; gains in 8p occurred only in conjunction with gains in 8q. The smallest regions of copy number alteration were mapped to gain of 8q21 and loss of 6q21, 10q21, and 16q22. Sublocalization of these chromosomal changes to single-band resolution should accelerate the identification of genes involved in the genesis of ocular melanoma.

16q24.1 microdeletion in a premature newborn: Usefulness of array-based comparative genomic hybridization in persistent pulmonary hypertension of the newborn.

OBJECTIVE:: Report of a 16q24.1 deletion in a premature newborn, demonstrating the usefulness of array-based comparative genomic hybridization in persistent pulmonary hypertension of the newborn and multiple congenital malformations. DESIGN:: Descriptive case report. SETTING:: Genetic department and neonatal intensive care unit of a tertiary care children's hospital. INTERVENTIONS:: None. PATIENT:: We report the case of a preterm male infant, born at 26 wks of gestation. A cardiac malformation and bilateral hydronephrosis were diagnosed at 19 wks of gestation. Karyotype analysis was normal, and a 22q11.2 microdeletion was excluded by fluorescence in situ hybridization analysis. A cesarean section was performed due to fetal distress. The patient developed persistent pulmonary hypertension unresponsive to mechanical ventilation and nitric oxide treatment and expired at 16 hrs of life. MEASUREMENTS AND MAIN RESULTS:: An autopsy revealed partial atrioventricular canal malformation and showed bilateral dilation of the renal pelvocaliceal system with bilateral ureteral stenosis and annular pancreas. Array-based comparative genomic hybridization analysis (Agilent oligoNT 44K, Agilent Technologies, Santa Clara, CA) showed an interstitial microdeletion encompassing the forkhead box gene cluster in 16q24.1. Review of the pulmonary microscopic examination showed the characteristic features of alveolar capillary dysplasia with misalignment of pulmonary veins. Some features were less prominent due to the gestational age. CONCLUSIONS:: Our review of the literature shows that alveolar capillary dysplasia with misalignment of pulmonary veins is rare but probably underreported. Prematurity is not a usual presentation, and histologic features are difficult to interpret. In our case, array-based comparative genomic hybridization revealed a 16q24.1 deletion, leading to the final diagnosis of alveolar capillary dysplasia with misalignment of pulmonary veins. It emphasizes the usefulness of array-based comparative genomic hybridization analysis as a diagnostic tool with implications for both prognosis and management decisions in newborns with refractory persistent pulmonary hypertension and multiple congenital malformations.

Microarray based real-time analysis of nucleic acid hybridization kinetics and thermodynamics

Ziel dieser Dissertation ist die experimentelle Charakterisierung und quantitative Beschreibung der Hybridisierung von komplementären Nukleinsäuresträngen mit oberflächengebundenen Fängermolekülen für die Entwicklung von integrierten Biosensoren. Im Gegensatz zu lösungsbasierten Verfahren ist mit Microarray Substraten die Untersuchung vieler Nukleinsäurekombinationen parallel möglich. Als biologisch relevantes Evaluierungssystem wurde das in Eukaryoten universell exprimierte Actin Gen aus unterschiedlichen Pflanzenspezies verwendet. Dieses Testsystem ermöglicht es, nahe verwandte Pflanzenarten auf Grund von geringen Unterschieden in der Gen-Sequenz (SNPs) zu charakterisieren. Aufbauend auf dieses gut studierte Modell eines House-Keeping Genes wurde ein umfassendes Microarray System, bestehend aus kurzen und langen Oligonukleotiden (mit eingebauten LNA-Molekülen), cDNAs sowie DNA und RNA Targets realisiert. Damit konnte ein für online Messung optimiertes Testsystem mit hohen Signalstärken entwickelt werden. Basierend auf den Ergebnissen wurde der gesamte Signalpfad von Nukleinsärekonzentration bis zum digitalen Wert modelliert. Die aus der Entwicklung und den Experimenten gewonnen Erkenntnisse über die Kinetik und Thermodynamik von Hybridisierung sind in drei Publikationen zusammengefasst die das Rückgrat dieser Dissertation bilden. Die erste Publikation beschreibt die Verbesserung der Reproduzierbarkeit und Spezifizität von Microarray Ergebnissen durch online Messung von Kinetik und Thermodynamik gegenüber endpunktbasierten Messungen mit Standard Microarrays. Für die Auswertung der riesigen Datenmengen wurden zwei Algorithmen entwickelt, eine reaktionskinetische Modellierung der Isothermen und ein auf der Fermi-Dirac Statistik beruhende Beschreibung des Schmelzüberganges. Diese Algorithmen werden in der zweiten Publikation beschrieben. Durch die Realisierung von gleichen Sequenzen in den chemisch unterschiedlichen Nukleinsäuren (DNA, RNA und LNA) ist es möglich, definierte Unterschiede in der Konformation des Riboserings und der C5-Methylgruppe der Pyrimidine zu untersuchen. Die kompetitive Wechselwirkung dieser unterschiedlichen Nukleinsäuren gleicher Sequenz und die Auswirkungen auf Kinetik und Thermodynamik ist das Thema der dritten Publikation. Neben der molekularbiologischen und technologischen Entwicklung im Bereich der Sensorik von Hybridisierungsreaktionen oberflächengebundener Nukleinsäuremolekülen, der automatisierten Auswertung und Modellierung der anfallenden Datenmengen und der damit verbundenen besseren quantitativen Beschreibung von Kinetik und Thermodynamik dieser Reaktionen tragen die Ergebnisse zum besseren Verständnis der physikalisch-chemischen Struktur des elementarsten biologischen Moleküls und seiner nach wie vor nicht vollständig verstandenen Spezifizität bei.

Microarray-based comparative genomic hybridisation of breast cancer patients receiving neoadjuvant chemotherapy

We analysed the molecular genetic profiles of breast cancer samples before and after neoadjuvant chemotherapy with combination doxorubicin and cyclophosphamide (AC). DNA was obtained from microdissected frozen breast core biopsies from 44 patients before chemotherapy. Additional samples were obtained before the second course of chemotherapy (D21) and after the completion of the treatment (surgical specimens) in 17 and 21 patients, respectively. Microarray-based comparative genome hybridisation was performed using a platform containing approx5800 bacterial artificial chromosome clones (genome-wide resolution: 0.9 Mb). Analysis of the 44 pretreatment biopsies revealed that losses of 4p, 4q, 5q, 12q13.11–12q13.12, 17p11.2 and 17q11.2; and gains of 1p, 2p, 7q, 9p, 11q, 19p and 19q were significantly associated with oestrogen receptor negativity. 16q21–q22.1 losses were associated with lobular and 8q24 gains with ductal types. Losses of 5q33.3–q4 and 18p11.31 and gains of 6p25.1–p25.2 and Xp11.4 were associated with HER2 amplification. No correlations between DNA copy number changes and clinical response to AC were found. Microarray-based comparative genome hybridisation analysis of matched pretreatment and D21 biopsies failed to identify statistically significant differences, whereas a comparison between matched pretreatment and surgical samples revealed a statistically significant acquired copy number gain on 11p15.2–11p15.5. The modest chemotherapy-driven genomic changes, despite profound loss of cell numbers, suggest that there is little therapeutic selection of resistant non-modal cell lineages.

Transcriptional Response of Peripheral Lymphocytes to Early Fibrosarcoma: A Model System for Cancer Detection Based on Hybridization Signatures

Since circulating leukocytes, mainly B and T cells, continuously maintain vigilant and comprehensive immune surveillance, these cells could be used as reporters for signs of infection or other pathologies, including cancer. Activated lymphocyte clones trigger a sensitive transcriptional response, which could be identified by gene expression profiling. To assess this hypothesis, we conducted microarray analysis of the gene expression profile of lymphocytes isolated from immunocompetent BALB/c mice subcutaneously injected with different numbers of tumorigenic B61 fibrosarcoma cells. Flow cytometry demonstrated that the number of circulating T (CD3(+)CD4(+) or CD3(+)CD8(+)) or B (CD19(+)) cells did not change. However, the lymphocytes isolated from tumor cell-injected animals expressed a unique transcriptional profile that was identifiable before the development of a palpable tumor mass. This finding demonstrates that the transcriptional response appears before alterations in the main lymphocyte subsets and that the gene expression profile of peripheral lymphocytes can serve as a sensitive and accurate method for the early detection of cancer. Exp Biol Med 234:802-812, 2009

Microarray-based detection of 90 antibiotic resistance genes of gram-positive bacteria

A disposable microarray was developed for detection of up to 90 antibiotic resistance genes in gram-positive bacteria by hybridization. Each antibiotic resistance gene is represented by two specific oligonucleotides chosen from consensus sequences of gene families, except for nine genes for which only one specific oligonucleotide could be developed. A total of 137 oligonucleotides (26 to 33 nucleotides in length with similar physicochemical parameters) were spotted onto the microarray. The microarrays (ArrayTubes) were hybridized with 36 strains carrying specific antibiotic resistance genes that allowed testing of the sensitivity and specificity of 125 oligonucleotides. Among these were well-characterized multidrug-resistant strains of Enterococcus faecalis, Enterococcus faecium, and Lactococcus lactis and an avirulent strain of Bacillus anthracis harboring the broad-host-range resistance plasmid pRE25. Analysis of two multidrug-resistant field strains allowed the detection of 12 different antibiotic resistance genes in a Staphylococcus haemolyticus strain isolated from mastitis milk and 6 resistance genes in a Clostridium perfringens strain isolated from a calf. In both cases, the microarray genotyping corresponded to the phenotype of the strains. The ArrayTube platform presents the advantage of rapidly screening bacteria for the presence of antibiotic resistance genes known in gram-positive bacteria. This technology has a large potential for applications in basic research, food safety, and surveillance programs for antimicrobial resistance.

Parallel human genome analysis: microarray-based expression monitoring of 1000 genes.

Microarrays containing 1046 human cDNAs of unknown sequence were printed on glass with high-speed robotics. These 1.0-cm2 DNA "chips" were used to quantitatively monitor differential expression of the cognate human genes using a highly sensitive two-color hybridization assay. Array elements that displayed differential expression patterns under given experimental conditions were characterized by sequencing. The identification of known and novel heat shock and phorbol ester-regulated genes in human T cells demonstrates the sensitivity of the assay. Parallel gene analysis with microarrays provides a rapid and efficient method for large-scale human gene discovery.

A study of adrenocortical tumors in children and adolescents by a comparative genomic hybridization technique

Adrenocortical tumors (ACT) are rare neoplasms of the adrenal glands accounting for 0.2% of all pediatric cancers. However, the incidence of ACT in South Brazilian children is 10 to 15 times greater than the worldwide incidence. Comparative genomic hybridization studies have revealed the presence of a high degree of chromosomal instability in ACT. We evaluated 16 ACT, 8 of them carcinomas and 8 adenomas. The presence of changes in DNA copy numbers was determined by comparative genomic hybridization, and the findings were validated by real-time polymerase chain reaction on the basis of IGF-II gene expression. The adenomas showed a mean of 19.7 imbalances per case, with the most frequent gain and loss being 4p15.1-p15.3 and 20p11.2-p13.2, respectively. The carcinomas presented with a mean of 35.5 imbalances per case, with the more frequent gain being 2q14.1-q24.3 and the more frequent losses being 3q21-q26.2, 20q12-qter, and 22q11.2-q13.3. The most frequent imbalances in both adenomas and carcinomas were gains of 1p21-p31.2, 2p12-p21 and loss of 20p11.2-p12. The expression of IGF-II mRNA (11p15.5) was higher in samples that presented with a gain of this region. It has been established that great genomic instability exists in pediatric ACT.

Reliability of short comparative genomic hybridization in fibroblasts and blastomeres for a comprehensive aneuploidy screening: first clinical application

BACKGROUND: Comparative genomic hybridization (CGH) is a valuable alternative to fluorescence in situ hybridization (FISH) for preimplantation genetic screening (PGS) because it allows full karyotype analysis. However, this approach requires the cryopreservation of biopsied embryos until results are available. The aim of this study is to reduce the hybridization period of CGH, in order to make this short-CGH technique suitable for PGS of Day-3 embryos, avoiding the cryopreservation step. METHODS: Thirty-two fibroblasts from six aneuploid cell lines (Coriell) and 48 blastomeres from 10 Day-4 embryos, discarded after PGS by FISH with 9 probes (9-chr-FISH), were analysed by short-CGH. A reanalysis by the standard 72 h-CGH and FISH using telomeric probes was performed when no concordant results between short-CGH and FISH diagnosis were observed. The short-CGH was subsequently applied in a clinical case of advanced maternal age. RESULTS: In 100% of the fibroblasts analysed, the characteristic aneuploidies of each cell line were detected by short-CGH. The results of the 48 blastomeres screened by short-CGH were supported by both 72 h-CGH results and FISH reanalysis. The chromosomes most frequently involved in aneuploidy were 22 and 16, but aneuploidies for the other chromosomes, excepting 1, 10 and 13, were also detected. Forty-one of the 94 aneuploid events observed (43.6%) corresponded to chromosomes which are not analysed by 9-chr-FISH. CONCLUSIONS: We have performed a preliminary validation of the short-CGH technique, including one clinical case, suggesting this approach may be applied to Day-3 aneuploidy analysis, thereby avoiding embryo cryopreservation and perhaps helping to improve implantation rate after PGS.

Oral squamous cell carcinoma: status of tight junction claudins in the different histopathological patterns and relationship with clinical parameters. A tissue-microarray-based study of 136 cases

Aims Claudins are integral transmembrane proteins of the tight junctions, critical for maintaining cell adhesion and polarity. Alterations in the expression of individual claudins have been detected in carcinomas and appear to correlate with tumour progression. Methods In this study, a panel of anti-claudin antibodies (anti-claudins 1, 2, 3, 4, 5 and 7) was employed to map claudin expression in 136 cases of oral squamous cell carcinoma (OSCC) organised in a tissue microarray. Results Claudins were expressed in a reticular pattern up to the prickle layer in normal mucosal epithelium. In OSCC, claudins were strongly present in well-differentiated tumours, they presented mild and low expression in moderately differentiated OSCC, and were negative in poorly differentiated OSCC; the absences of claudin 1 (p = 0.002) and claudin 4 (p<0.001) were associated with moderately/poorly differentiated tumours. Strong expression of claudin 4 was associated with decreased perineural infiltration (p = 0.024). Claudins 5 and 7 were mostly negative or weakly expressed in all cases studied. Expression of claudin 7 was associated with the early clinical stages of the disease, whereas loss of claudin 7 tended to be more frequent in advanced stages of OSCC (p = 0.054). Absence of claudin 7 was also associated with absent vascular infiltration (p = 0.045) and with presence of recurrence (p = 0.052). Conclusions Claudin expression patterns showed a strong correlation with histological type of OSCC; claudin expression was decreased in areas of invasion, and negative in poorly differentiated tumours. This pattern may be related to evolution and prognosis of these tumours, especially in the case of claudin 7, which seems to be associated with a poor prognosis in OSCC.

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

A microarray-based system for the simultaneous analysis of single nucleotide polymorphisms in human genes involved in the metabolism of anti-malarial drugs

Background: In order to provide a cost-effective tool to analyse pharmacogenetic markers in malaria treatment, DNA microarray technology was compared with sequencing of polymerase chain reaction (PCR) fragments to detect single nucleotide polymorphisms (SNPs) in a larger number of samples. Methods: The microarray was developed to affordably generate SNP data of genes encoding the human cytochrome P450 enzyme family (CYP) and N-acetyltransferase-2 (NAT2) involved in antimalarial drug metabolisms and with known polymorphisms, i.e. CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, and NAT2. Results: For some SNPs, i.e. CYP2A6*2, CYP2B6*5, CYP2C8*3, CYP2C9*3/*5, CYP2C19*3, CYP2D6*4 and NAT2*6/*7/*14, agreement between both techniques ranged from substantial to almost perfect (kappa index between 0.61 and 1.00), whilst for other SNPs a large variability from slight to substantial agreement (kappa index between 0.39 and 1.00) was found, e. g. CYP2D6*17 (2850C>T), CYP3A4*1B and CYP3A5*3. Conclusion: The major limit of the microarray technology for this purpose was lack of robustness and with a large number of missing data or with incorrect specificity.